1 



April 30, 1998 UMN-TH-1610/98 



00 ■ 

o^ ■ Bosonized Massive N-flavor Schwinger Model 

u 

Oh: 
< 
o _ 

^ Yutaka Hosotani^ and Ramon Rodriguez^'^ 

1^ I ^School of Physics and Astronomy, University of Minnesota 

O ■ Minneapolis, MN 55455-0112, USA 

(N : 

Tt . ^SDRC, Minneapolis, MN 55425-1528, USA 

O 

00 



Oh! 

r-| ! Abstract 



X 



The massive N-flavor Schwinger model is analyzed by the bosonization 
method. The problem is reduced to the quantum mechanics of N degrees 
5^ I of freedom in which the potential needs to be self-consistently determined by 

its ground-state wave function and spectrum with given values of the 9 pa- 
rameter, fermion masses, and temperature. Boson masses and fermion chiral 
condensates are evaluated. In the N=l model the anomalous behavior is found 
at ^ ~ vr and m/^u ~ 0.4. In the N=3 model an asymmetry in fermion masses 
(mi < m2 ^ m^) removes the singularity at 9 = tt and T = 0. The chiral 
condensates at 9 = are insensitive to the asymmetry in fermion masses, but 
are significantly sensitive at 9 = n. The resultant picture is similar to that 
obtained in QCD by the chiral Lagrangian method. 
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1. Introduction 

Two-dimensional quantum electrodynamics (QED2) "witli massive N-fiavor 

fermions resembles with four-dimensional quantum cliromodynamics (QCD) in many 
respects. In both theories non- vanishing chiral condensates are dynamically generated. 
Fractionally charged test particles are confined in QED2, whereas quarks, or colored objects, 
are confined in QCD. The dynamical chiral symmetry breaking and confinement are not 
independent phenomena in QED2, however. There would be no confinement if there were 
no chiral condensates. 0-0 

It has also been recognized that QED2 describes spin systems in nature. § A spin ^ anti- 
ferromagnetic spin chain is equivalent to two-flavor massless Schwinger model in a uniform 
background charge density. A n-leg spin ladder systems is equivalent to a coupled set of n 
Schwinger models. This equivalence has been successfully employed to account for the gap 
generation in spin ladder systems. 0] 

QED2 has been analyzed by various methods. On the analytic side it has been inves- 
tigated in the perturbation theory, in the path integral method, and in the bosonization 
method. In a series of papers we have shown how to evaluate chiral condensates and boson 
mass spectrum for arbitrary values of the 6 parameter, fermion masses (m), and temperature 
(T) by bosonization. |T^, [T^, |T2[ The mass perturbation theory has been formulated in the 
one-flavor model. [Il3|, [14| 

Investigation in the light-cone quantization method has been pushed forward both on 
the analytic and numerical sides. |15|]-[0 The bound state spectrum has been evaluated in 



the entire range of a fermion mass at 6' = and T = 0. Subtleties in the chiral condensate 
in this formalism has been noted. 

There has emerged a renewed interest in QED2 in the lattice gauge theory approach 
as well.|^T|-||23| Recently extensive simulations have been carried out for the = 1 and 



N = 2 models. Chiral condensates in the A^ = 1 model were evaluated at = and T = 
up to m/e < 1. The boson mass in the N = 2 model was evaluated for m/e < .5. After 
subtracting condensates in free theory (e = 0), which depends on regularization methods 
employed, one flnds a modest agreement between the bosonization and lattice results. 



In this paper we further exploit the bosonization method to investigate the dependence 
of chiral condensates and boson mass spectrum on the 6 parameter, fermion masses, and 
temperature. The advantage of our method lies in the ability of evaluating physical quantities 
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for arbitrary values of the 6 parameter and temperature. The current method, however, 
involves an approximation which is not valid for large fermion masses. Improvement is 
necessary in this direction. 

The Lagrangian of the model is given by 

C = ~\ F,,F^^ + Y: [i^a^iid, - CaA,)^, - mjy^l - m:^^2V^^} (1.1) 

a 

where '?/'2 = |(1 — 'y^)4'a and ipR = |(1 + 7^)'?Aa- Each field carries a charge and mass m^. 
We analyze the model on a circle (S^) with a circumference L. The boundary conditions are 
specified by 

A^{t,x + L) = Af,{t,x) 

^„(t,x + L) = -e2™'"^„(t,x) . (1.2) 

It is important to recognize that from the analysis on 5*^ one can extract physics at 
finite temperature T defined on a line R^. In the Matsubara formalism of finite temperature 
field theory, boson and fermion fields obey periodic or anti-periodic boundary condition in 
the imaginary time axis, respectively; 

A^^{T + ^,x) = A^{t,x) 

1pa{T + ^,X) = -iJa{T,x) . (1.3) 

Hence, if one, in a theory defined on 5^ with the boundary conditions eta = in ( |1.2| ), 
analytically continues t from the real axis to the imaginary axis, and then interchanges (or 
relabels) it and x, one arrives at a theory which is exactly the same as the theory defined 
on i?^ at T = L~^. This is a powerful equivalence. One can evaluate chiral condensates, 
Polyakov loops, and various correlators at T 7^ with the aid of this correspondence. 

This paper is organized as follows. In Section 2 the bosonized Hamiltonian is derived 
on S^. In Section 3 the 6 vacuum is introduced and the equation satisfied by its wave func- 
tion is derived. In Section 4 we show how the boson mass spectrum and chiral condensates 
are evaluated. Sections 3 and 4 together form a basis of our formulation. It defines the 
generalized Hartree-Fock approximation. The rest of the paper is devoted to applying the 
generalized Hartree-Fock equation to various models to evaluate the boson spectrum and 
chiral condensates. The case of massless fermions is discussed in Section 5. A useful trun- 
cated formula is derived in Section 6. The detailed analysis of the massive one-flavor model 
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is given in Section 7. The multi-flavor model with degenerate fermion masses is analyzed in 
Section 8. The case of general masses in the N=2 and N=3 models is investigated in Section 
9. A brief summary is given in Section 10. Three appendices collect useful formulas. 



2. The bosonized Hamiltonian 

Our basic tool is the bosonization method, with which we shall reduce the model (|1 . 1| ) 
to a quantum mechanical system of flnite degrees of freedom. The bosonization method has 



been developed in many body theory, 1 24] and in field theories on a line |]25| , 0, 26]. The 
method has been elaborated on a manifold in the context of string theory. |27|] 

The bosonization on S"^ is particularly unamabiguous, sustaining the absolute rigor. 
QED2 on was first studied by Nakawaki PB[ and has been developed by many authors 



PT| ] It was simplified in ref. [^, which we follow in the sequel. 



In this section we present a brief review of the method, applying it to the system (|1.1|) . 
Although the essence is well known in the literature, subtle factors associated with the 
multi-fiavor nature and implementation of arbitrary boundary conditions are worth spelling 
out. With clever choice of Klein factors the Hamiltonian of the iV-fiavor Schwinger model is 
transformed into a surprisingly simple bosonized form. 



We note that the A^-fiavor Schwinger model has been analysed in refs. p6|, |10|, |Il|, p^j . 
The model at finite temperature, which is equivalent to the model on S^, has been analysed 
in ref. ^3[. The model on other manifolds also have been investigated in ref. |3^. The 



conformal field theory approach to QED2 has been proposed by Itoi and Mukaida [p5|] , 
which has many features in common with our bosonization approach. Correlators of various 



physical quantities have been discussed in ref. . 

Bosonization of an arbitrary number of fermions on a circle (0 < x < L) obeying 
boundary conditions 

V'a(t,a: + L) = -e2™'^Va(t,a;) {a = 1 ^ N) (2.1) 

is first carried out in the interaction picture defined by free massless fermions: i'ydip = 
0.|2^, ^ We introduce bosonic variables: 

ab ra b,] 
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all other commutators = 



OO -| 



In terms of these variables = ('?/'!^, can be expressed as [P7 



(2.2) 



(2.3) 



where the Klein factors are given by |]T0| 

CI = exp|^7rX:(pt + /--2a6)| 
b=i ^ 

^ b=i ^ 



(2.4) 



CI 



This choice of Klein factors turns out very convenient. They satisfy 

ip±{x) = sign (6 > a) ^±{x) 

ip±ix) = sign {b > a) ^Ij^{x) 



C 



■at 



CI 

Cl^ 



(2.5) 



where sign (A) is defined to be +1 (—1), when A is true (false). By construction [dt =F 
dr^ ip^it, x) = 0, i.e. ipa satisfies a free massless Dirac equation: {•y^do + •y^di) tp"- = where 
Y = {cri,i(^2) and 7^ = 7%^ = -0-3. 

Under a translation along the circle 



^^(t, x + L) = -e'^'P± ^l{t, x) = -ij^it, x) e'^'Pi 



so that the boundary condition ( p..2|) is ensured by a physical state condition 



(2.6) 



e^'^Pi I phys ) = e^"^"" | phys ) . 



(2.7) 



Further conditions can be consistently imposed on physical states such that the Klein factors 
act in physical space as a semi-identity operator: | phys) = (+ or — )| phys). We shall 
see below that the Hamiltonian commutes with p°, — . 
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The fields {ip±{x)} satisfies desired equal-time anti-commutation relations. With the 
aid of ( pTsD , and ( p.2| ), it is straightforward to show 



others = (2.8) 

where a,f3 = + or — . Notice that the extra phase factors in ( p.8|) manifest the translation 
property ( p.6|) . The bosonization in the interaction picture is defined by ( p.3|) and ( p.TI) . 



In applying the bosonization method to the model (|1 . 1|) , it is most convenient to take 
the Coulomb gauge 

A^{t,x)=b{t) , 

Ao{t, x) = - dy G{x - y) jly^{t, y) , jg^ = H e^V'I^a , 

G{x + L)=G{x) , ^G{x) = 5l{x)-^ (2.9) 

in which the zero mode h{t) of Ai{t, x) is the only physical degree of freedom associated with 
the gauge fields. The Hamiltonian is given by 

1 /"^ r— — — 1 

^ ^ 2Z^b + dxY^ [i^al^-idl + eah)ll)a + mJjaLi^aR + ml^jj aR^aLj 

1 

-- / dxdy j^^{x)G{x - y)j^M{y) . (2.10) 
2 JO 

Here Pb is the momentum conjugate to b: Pb = Lb. The anti-symmetrization of fermion 
operators is understood. 

At all stages of bosonization, the gauge invariance must be maintained. Due caution is 
necessary in bosonizing a product of two field operators at the same point as it has to be 
regularized. At equal time the bosonization formula (p^.3D leads to 



Gauge invariant regularization amounts to dropping the P/{x — y) term. 



7 



Hence we have 



= - 1^ + i ^ + ^ ^^^^^^ ^^^^^^ 

In particular, currents are given by 

jl = -wrlr^] + - + e„- + ^a,^^ (2.13) 

where 0a = + (j)"L. In terms of 0^ = 4>a+ + 0a- where (pa± = (47r)~^/^ [g^ + 27rp^ (t ± 
a;)/L] + 0^, the current takes a simpler form 

j^^-^e^''d.4>a + S^'^ . (2.14) 

In the following discussions, however, we shall find that treating the zero mode and oscillatory 
mode parts separately is more convenient. 

The kinetic energy term is transformed to 

= + ^ : + e„6 + V^dAXf + - Vl^d^r-f : ■ (2.15) 

Putting all things together, we find 

H — Hq + + -ffmass 







a=l 
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■'^ ^ a=l 

/^^ = -Ee^-- . <^ = E-<l>a (2.16) 

where -ffmass represents the fermion mass term. It is convenient to express the Hamiltonian 
in terms of 

qa = ql + q- , Pa = |(p++p") 
^a = |(g+-g") , Pa = Pl-P- 

[qa,Pb\ = [qa,Pb\ = iSab , all others = ; (2.17) 

Hq becomes 

AT p2 1 

A few important conclusions can be drawn in the massless fermion case (-ffmass = 0). 
The zero mode part Hq and oscillatory part decouple from each other. Each part is 
bilinear in operators so that the Hamiltonian is solvable. The oscillatory part consists of one 
massive boson (0) with a mass n and — 1 massless bosons. The zero mode part must be 
solved with the physical state condition 

g^^lphys ) = - E e^plphys ) = . (2.19) 

a 

3. 9 vacuum 

When all ratios of various charges are rational, there results a 6 vacuum. In this 
article we restrict ourselves to the case in which all fermions have the same charges: = e. 
It is appropriate to introduce the Wilson line phase Qw- 

Qw = ebL , e*0^(*) =e''=/o^'^-^i(*'^) . (3.1) 
The zero mode part of the Hamiltonian (|2.18|) becomes 

^^0 = + + ^ E {{Qw + 2np.r + ^pI) (3.2) 
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where /i^ — Nc^/tt and Pw — Qw/^^L is the conjugate momentum to Qw In terms of the 
new variables 

, 27r ^ ^ , 27r 

[Q'^r,Pw]=i , [q'a,Pb\ = iSab , others = (3.3) 
the Hamiltonian becomes 

There appears an additional symmetry {QwiPa) ~^ {~^w-i ~Pa) when = and |. The 
Hamiltonian is invariant under 

©ly ^ ©w^ + 27r , Pa ^ Pa - 1 (3.5) 

or equivalently, in terms of the original fields, 

1 27rT 

The transformation is generated by a unitary operator 

[U,H] = . (3.7) 

In a vector-like theory Pa= P^— P- takes integer eigenvalues. Further [pa, H] = 0. We 
can restrict ourselves to states with p^ = as the energy is minimized there. The vacuum 
state is written as a direct product of ground states of the zero mode sector and oscillatory 
(0) mode sector. The ground state in the oscillatory mode sector is defined with respect 
to physical boson masses /laS. As we shall see, the ground state wave function in the zero 
mode sector affects the physical boson masses, and vice versa. These two must be determined 
self-consistently. Note that if there is a background charge {Qb.g. 7^ 0) as in the case of spin 
chains, then -Y^^aPa = -Qb.g.- 

With this understanding the vacuum wave function is written as 

/oo _ 
dpw X] \Pw,n,ra) f{pw,n,ra) (3.8) 
-00 f 1 

{n,ra} 



10 



Bosonized Schwinger model 



where \pw, n, ) is an eigenstate of Pw and pa- 
Pw\Pw,n,ra) = Pw\Pw,n,ra) 



Pa\Pw,n,ra) - i + hTa = N. 



It follows that 

e'''^'^ \pw,n,n) = \pw + k,n,n) 



e^- \p^, n,n) = j !^^' ^' ± ^"'l^^ , ';'-<l (3.10) 

1 \pw,n ± l,rbT i^} tor a = A/. 



Since U in (|3.7|) commute with the Hamiltonian, one can take 

f/|^vacW)=e+^^|v&.acW) • (3.11) 

As U\pw,n,ra) = e^'''P^\pw,n + l,ra) , one finds f{pw,n,ra) = e-^"''+2^^("+'^W /(^^^ 
where a = Yla=i ^a- 

This is the Q vacuum. The existence of the operator U was suspected long time ago 
in ref. 0, though its explicit form was not given. On a circle V is unambiguously written 
in terms of Piy aiid g^'s- This definition was first given in ref. |3D|. In passing, it has been 



noticed recently that similar definition of the Q vacuum arises in the framework of light cone 



quantization of the Schwinger model. |18, 19 



It is convenient to adopt a coherent state basis given by 

n,^a)= (^l,y, E e-('^^^^+-+'---^--) \pw, n,r^) . (3.12) 



Then 



fiPw.^a) = .(Li)/2 Eg^^'°'^''/(PH/,rJ . (3.13) 

^ {ra} 

The normalization is ( ^vac(^') l^vac(^) ) = 52^(0' - 9) and J dpw[dip] |/T = 1- 

The function f{pw, fa) is determined by solving the eigenvalue equation iftotl^vac(^) ) = 
-E'vacl^vac(^) )• Wc writc, for an operator Q = Q{Qw, Pw-,Pa)-, 
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Noticing Qw\pw-,n^'^a) = —i{d/dpw)\pw,n,i^a), one finds that 

Q = Qw + -^T.Pa ^ Q = 4j^--^T.^\ ■ (3-15) 

N ^1 I dpw N ^1 dfa J 

The operator Hq corresponding to Hq is 

Here the 99-Laplacian is given by 

A. ^ E (5^ - .^.)^ - ^ E - .^.) (I; - ./S.) (3.1T) 

where /3a = — «Af- 

The mass operator in the Schrodinger picture is 

Mf„ = = -Cl^Cl ■ 6-2-*^"^/-^ e^'^" A/'o[e*^'^"] (3.18) 

where the normal ordering No[ ] is defined with respect to massless fields. In the presence 
of -f^rnass, all boson fields (0) become massive. We denote a mass eigenstate with a mass fia 

by Xa- 

Xa = U^a(l>a , U'U = I . (3.19) 

The vacuum is defined with respect to these Xa fields. 
With the aid of ( [B.7| ), the mass operator becomes 

N 

n 

0=1 

N 



0=1 

5, = f{B{^^^L)^''-'^^' (3.20) 



Further 



so that 



Af-i f 5n> n e^''^" for a < AT 

= 5(pW - Pw^) n '^2.(</^; - </^6) ' (3.21) 
^=1 Un',n±i e^*^""" fora = A^ 

giiv^a for a < iV 

g = e^^'?'^ ^ Q = <( (3.22) 



12 Bosonized Schwinger model 

Let us write a fermion mass as nia — \ma \ e*'^" and drop the absolute value sign hence- 
forth. Then i^mass — ! dx ma(Mae*'^" + h.c. ). When i^mass acts on|^'vac(^) ), in general 
Xa quanta are also excited as well as in the zero-mode sector. In deriving an equation for 

f{pw,'^), we ignore those x-excitations, with the understanding that physical masses /Xq's 
are taken. Then 

Af-l 

-^mass = - XI 2maBa COs{ipa + So) - 2mNBN COs{9 -"^(fib- ^TTpw + Sn) ■ (3.23) 

0=1 

f{PW: <Pa) must satisfy HtotfiPW, <Pa) = Ey^fiPW: <Pa) whcrC ^tot = -^0 + -^mass- 

At this stage we recognize that it is appropriate to introduce 

f{PW, ^a) = KpW, - - Sa) . (3.24) 



The eigenvalue equation Htotfipw, '-Pa) = E^acf{pwi '-Pa) now reads 

|-(2^j ^-(A^-l)A^ + nPH^,(^)|/(pty,(^) = e/(p^^,(^) . (3.25) 
where e = {NLE^/2'k) + (-kN'^/12). The potential is given by 

V{PW: <P) = H ^ Pw maLBa COS [(pa ^ j 

^ a=l 

N-1 N 
^ O^S - ^ iPa , OeS = + J2^a ■ (3-26) 
a=l a=l 

The vacuum is 

|*vac(^) ) = I dpw[d^p\ \PW, n, ipa ) 

v27r n •> 

X e— f{pw. Va + ^ + K) . (3.27) 
The problem has been reduced to solving the Schrodinger equation for N degrees of freedom. 
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4. Boson masses and condensates 



In deriving the equation ( |3.25| ) for the vacuum wave function f{pw,ip), we have sup- 



posed that we already know the masses /ia's of the boson fields Xa's- In this section we show 
how these related to the vacuum wave function f{pw, f) itself. Hence we obtain a 

self-consistency condition for the vacuum. Along the way we shall also find that the chiral 
condensate {ipai^a) is related to /ia's. 
From (|3.22|) it follows that 

(e±*^'"), = lim(vl/.,e(^')|e^'''1*vacW)/(^vac(^')l^vacW) 



where v^at is defined in ( p.26| ) and the /-average is given by ( F )j = / dpw[dip] In our 

approximation scheme {Maa)e = —{Ba/L){e^'^"- )g, and therefore 

2— 27rpw 

-Ba{ COS{ipa — 



[-ipa-ipaYe = -T^o,{ COs{ipa ^ " ^a) )/ • (4-2) 



Rigorously speaking, the formula (|4.2| ) is valid only for small m. Adam has determined 
the condensate to 0(m) in mass perturbation theory in the = 1 case.[|T^ Although the 



formula ( [4.2|) incorporates some of the effects of a fermion mass m through Ba and f{(p), it 
does not incorporate higher order radiative corrections which become important for a large 
m ^ fi. We shall see that our formula gives a fairly good agreement with the lattice result 
for m < /i in the A^ = 1 case. Adam's mass perturbation theory fails for m > .SjH.pO 



There is ambiguity in the definition of the composite operator ipiplx). It diverges in 
perturbation theory. It has to be normalized such that ( ■?/') = in a free theory (e = 0) 
in the infinite volume limit {L —>■ oo). In other words 

{^a^a)e= a^a Ye - a^a )Le ■ (4-3) 

The values of {ipai'a Ye and {ip ai'a )free depend on the regularization method employed, 
but the difference does not. In our regularization scheme we shall find {ip a'^aYiree = 
-e'^mjTi. (See (p|).) 

The fermion mass term -ffmass = J dx J2a {^a^^^" Maa + (h.c.) } has many effects. In 
addition to giving a "potential" in the zero mode sector as discussed in the previous section, 
it also gives mass terms (oc x^) and other interactions. It follows from ( p.20|) that 

i^mass ^ -T / E ^aBai e^(^-2w)/iV )^ -Q [^.l/.^v^^x.] + (h.^.) . (4.4) 

^ a a=l ^ 
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When fermion masses are small compared with the coupling constant, it is legitimate 
to expand //mass in power series of Xa- We define 



Ra +lla= ^ rUaBa " { 6 



i{ipa-2Trpw /N) \ 

If ■ 



(4.5) 



It follows that 



Including the additional mass term coming from the Coulomb interaction, one finds 



mass 



K, 



al3 



^ UggUlSb + UggUffaRa ■ 
a,b a 



(4.6) 



(4.7) 



Uaa^ are determined such that K^fs = fJ'aSa/s- In other words, we diagonahze 



^ = iV 



Ri 



Rn I 



K 



UlCU' 



V 



(4.8) 



The set of equations, (|3.2CI|) , (|3.25|) , ( [4.5|) , and ( [4. 81) needs to be solved simultaneously. This 
is a Hartree-Fock approximation applied to the zero mode and oscillatory mode sectors. We 
call it the generalized Hartree-Fock approximation 

In terms of Ra the chiral condensates are, in case 6a = 0, ma{ip ai^a )'e = — -Ra/4vr. Eqs. 
( [4.8| ) and ( ^4.3| ) relate boson masses to chiral condensates. It is a part of the PCAC (Partially 
Conserved Axial Currents) relations. 

As fermion masses become larger, nonlinear terms (~ x^) in -f^mass become relevant. 
The boson masses are not simply given by (|4.^ ). Improvement is necessary. 
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5. Massless fermions (aa = a) 

When all fermions are massless and satisfy the same boundary conditions = a, the 
model is exactly solvable. In this case Eq. ( |3.25| ) reduces to 

N-l Q2 2 ^1 



A 



The ground state, or the vacuum, wave function is independent of ipa- It is given by 

/(p^, y,) = const ■ e-"^^^'^/'^ . (5.2) 

The boson mass spectrum is given by /^i = = ^/Ne/Ti and /i2 = ■ ■ ■ = /^at = 0. The 
orthogonal matrix U in (|3.19|) has Uia = l/"\/iV and Ba in ( p.20| ) is 

B^ = B = B{^iLf"' . (5.3) 

For N > 2, { cos{ipa — 2TTpw/N))jr = in (|4.2| ) as f{pw,f) is independent of ipa- The 
chiral condensate {ipa'4^a)e vanishes for N > 2, reflecting Coleman's theorem which states 
that in two dimensions a continuous symmetry cannot be spontaneously broken. ||37[1 The 
non- vanishing {ipa'ipa)e breaks the SU{N) chiral symmetry. For = 1 the U{1) chiral 
symmetry is broken by anomaly. Recalling = 6' in ( |3.26|) for = 1, we have 

(*A)« = ( -I ^f"^' (5.4) 
[ for iV > 2. 

The condensate for = 1 is plotted in fig. 1. 

The SU (N) invariant condensate is non- vanishing, however. It follows from (|3.9|) that 

N N-l 

{P'w^ n', V?'| exp ( ± 2 ^ Qa) \pw, n,ip) = 6n',n±lS{Pw - Pw) n ^^A^'a " ^a) ■ (5.5) 
a=l a=l 

Straightforward calculations in the Schrodinger picture yield 

{^NL■■■^lL^l'lR■■■^l^NR)e = {^f {e-'^^''^No[e'^^'])g 

In the infinite volume limit it approaches e~*^( — iie'^ /Air)^ . 



e-^'e"^"/'^^ . (5.6) 



16 



Bosonized Schwinger model 




T/[i 



Figure 1: Temperature dependence of the chiral condensate in the A'^ = 1 model at ^ = 0. There 
appears crossover transition around T/fi = 1. 

6. Truncation 



There are two potential terms in the equation for f{pw,<^), Eq. (|3.25| ). If fermion 
masses are small compared to /i, one of them, (/iL)^pp^//4, dominates over the other term 
in V{pw,'^)- The condition is rriaBa <^ fJ^^L. For nia <C /i, it is satisfied when /iL ^ m//i 
{T/fi <^ fi/m). 

In such situation Pw acts as a fast variable, whereas (faS act as slow variables. The 
wave function can be approximated by f{pw,'^) = J2'^o'^s{pw)fs{'^) ~ u{pw)f{'^) where 
u{pw) = uoipw) is 



2tiJ dp^ 



Pw \ u{pw) 



w 



47r 



u{pw) 



u 



'2 /2N 



(6.1) 



Then cosine term in V]^{pw,'-p) is approximated by 

2'npw 



cos 



j dpw cos [ipa - 



l-np 



w 



N 



u{pwf 



(6.2) 



We remark that the truncated equation has symmetry py/ —pw, though the original 
equation does not in general. The truncation may not be good for physical quantities 
sensitive to this symmetry. 
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The function /(v?) satisfies 



-(iV-l)A^ + V5v(^)^/(^) = eo/(^) 
NL 



TT 



N 
a=l 



(6.3) 



This equation has been extensively studied in the N=2 and A^=3 cases in Refs ||I0|, |Tl|]. 
Let us denote )) j = J[dip] F{ip) |/(<^)p. Then (|aD and gj) become 



Ra + Ha 



— rUaBa e 



-n/NfiL// iipa 



(6.4) 



All formulas in Section 4 remain intact with these substitutions. In particular, the chiral 
condensate satisfies {ip a'^aYe = —Ra/^T^^a- 



7. = 1 massive case 

(1) Generalized Hartree-Fock approximation 

With one fermion there is no degree of freedom. We write mi = m,5i = 0, ai = a. 
Recall that fi = e/y^ and (fi = Ocr. The vacuum wave function is determined by 

+ —^Pw cos{2ttpw - ^cfr) }f{pw) = ^f{pw) ■ (7.1) 



When m fi, the boson mass /ii must satisfy 

/i? = /i' + ( cos(27rp^y - ^) )^ . (7.2) 

As m becomes larger, the formula ( |7.2| ) needs to be improved to incorporate nonlinear 
effects in the fermion mass m. In particular, /xi = 2m + 0(e^/m) for m ^ /x, as the boson 
is interpreted as a fermion-antifermion bound-state. 

In determining a physical chiral condensate, one needs to subtract a condensate in free 
theory as discussed in Section 4. The "free" limit corresponds to the limit m ^ /i. At the 
moment we do not have reliable formulas which relate /ii, m, and {^jj'tpYg for m ^ fi. The 
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best we can do is to extrapolate ( [4 .21) and ( [7.2|) for large m to determine the subtraction 
term within our approximation. 



It follows from ( [7.1|) that in the week coupling e/m <^ 1 and infinite volume L ^ oo 
limits ( cos(27rpi4/ — 9)) = 1. If the formula { \l.2\) is employed, one obtains /xi = 2e'm. 

- / e^T 

(^^)free = (7-3) 

TT 



Combined with ([4.2| ), it gives 
The chiral condensate is therefore given by 



>^), = -^^^ (cos(2w-^))^ + — m (7.4) 



which we expect to be a good approximation for m < fi. We stress that within our approxi- 
mation the subtraction term is given by ( |7.3| ). In an exact treatment the boson mass in the 
weak coupling limit should be given by ni = 2m. To achieve this, one has to improve both 
(|]2D and (|73) consistently. 

In the massless case (m = 0), {ipip)^~^ = —2L^^B(fiL)e~'"^^^ cos 6. In the infinite 
volume limit or zero temperature limit, it approaches —{fie'^/27T) cos6. 

It may be of interest to apply a perturbation theory |12[ to ( |7.1| ) when m fi. Write 
( [7. 1|) in the form 



{Ho + V)\^) = E\^) 
Hq = 27rfiL(a^a + ^) 

V = —Kcos{2'Kpw — d) , K. = 4:7imLB{fiiL) (7.5) 



where the annihilation and creation operators are 

1 / . 1 d 



In terms of number eigenstates \n), |\I' ) is found to be, to 0{m/fi), 

1 °° 1 

|^) = |0)-— — E-I^)(^l^|0) ■ (7.7) 

Since cos(27rpvF - &) = ie-'^A'^(^e"*^e*"'*V*°'' + h.c.) where a = {2tx/^LY''^, it follows 

(c 
/ 

cos{2Tipw - 0))f = ^ ^ -|(0|cos(27rpvK -^)|n) 



that ( cos(27rpv^ - 6) )f = e"^/'^-^ cos^ and 
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27r/iL n ■ n! V^y* 

K ( /•27r//tL p2 _ 1 /•27r/uL p-z _ 1 

e-^-/'^^ / t/^^ - + cos2e dz- 



2iTiiL 



Tfl 

— (1 - cos 26) for /iL > 1 . (7i 
/i 



The boson mass /ii = /x + is obtained from (|7.2|): 

4-7r 

5^ = — mS(/iL)e-"/^^cos^ + 0(mVAi2) 



~ me''' cos 6^ for yuL ^ 1 . (7.9) 

The chiral condensate is given by 

(^^)' = --{b{iiL) + B\iiL)5iiL\e-^'>'^ cosO - -B{iiL){cos{2Tipw -e))f 
L ^ ^ L ■' 

= -^/icos6' + ^(-3 + cos26')m for /xL > 1 . (7.10) 
27r 47r 



With (|7;E|) 



/icos6'H (l + cos26')m. (7.11) 

2n An 



This result differs from the result in the mass perturbation theory for the reason explained 
above and should not be taken seriously. Our formalism, however, allows us to estimate 
{ipi'){ree: therefore the physical condensate {ip4')0- The mass perturbation theory is 
valid only for small m/fi <^ 1, whereas our formalism allows the numerical evaluation of 
various physical quantities even for m/fir^ 1. We shall see below a good agreement between 
ours and the lattice gauge theory in the range m/^ < 1. See the subsection (3) below. 

(2) Mass perturbation theory 

Corrections can be evaluated in a power series in m/e. This analysis was carried out at 
zero temperature by Adam.[n| We present the analysis at finite temperature. 

We illustrate the computation for the chiral condensate. Recall ipip = 
-B{^L)L-\e'iK+ + e~'iKJ) where K± = e±*v^'^'"'e^^^<^'^'. We have suppressed ir- 
relevant factors. In the invariant perturbation theory 

{^^)^^^ = -tm\^^y fd'x Y: (T[e-''We^'"'(°)]T[ir,(x)ir,(0)]), . (7.12) 

^ ^ ) J a,b=± 
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Making use of (^) and the identity ( e^°<^'+'(^)e^^'^^ ) = e'^'^I'^'+H^),'/'' '(?/)] ^ one finds 
(^^)« = _^5(/iL)V2-/^^y'd2a;{(e^-*^(^)-l) + cos2^^(e-^-^^(^')-l)} 

The disconnected component has been subtracted. Deforming the t-integral to the imaginary 
axis, one finds 



TC 

r2iT 



X 



J dr j dx I (e+^(-'-;^^/2-) _ 1) + cos 20 ^^-E{r,x■,^,L/2.) _ ^) I 



= [F+{fiL) + F_{fiL) cos20]m 
1 °° 1 

E(r,x;z) = -e-^" + 2y — e"''"" cosnx , z;„ = (n^ + ^2)1/2 _ (7^4) 

The coefficient F^{^L) is finite, whereas diverges logarithmically near r = x = 

0. The divergence is due to the 0(m) correction in the free (e = 0) theory, which must be 
subtracted to define the physical chiral condensate as explained in Section 4. Hence in the 
mass perturbation theory 

T 

(^^)^ = -— /icos0 + m[F+(/iL) -F_?"" + F_(/iL)cos20l + •■■ . (7.15) 
27r L J 

F^'''^*^ is obtained in a similar manner. One starts with a massless free fermion theory. In 
the bosonization method the vacuum satisfies p|vac) = 0. Since ( e^*'' ) = 0, the condensate 
vanishes; (tptp) = 0. The fermion mass is treated as a perturbation. To 0(m) one gets an 
expression which is the same as ( |7.12| ) except that B{fiL) is replaced by 1 and (f){x) represents 
a massless field. Making use of (|C.2|), one finds that 



L2 

where Gfreei^) = ^i{T[(f>{x)(f){0)] ) is a massless propagator excluding the contribution from 
the zero mode. Employing ( |B.1| ) and deforming the integration path, one finds 



^ 777 /"cxj /'27r p 

{^^lj)\l>=-l!l dr dx- ^- —=mF^''. (7.17) 



Notice that F^'^'^ is independent of L. Comparing (|7.17|) with ( [7. 14|) , we observe that F_ 



free 

+ 
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F^{^L) can be easily evaluated by numerical integration. At fiL /2n = 10 ~ 50, -F_ = 
0.357. This is consistent with the number, 0.3581, obtained by Adam in the L —>■ oo limit. 
The divergence in F+{fiL) and F^^'^ makes the evaluation of the difference very difficult. 
We comment that Adam's subtraction procedure to get finite F^'^'^ is inconsistent. Indeed, 
his massless propagator differs from that obtained by taking the /i — limit of the massive 
propagator in ref. [|l^]. Further, his numerical estimate, —0.39126, for F-^.{oo)—F^^^ disagrees 
with the lattice result. At 6* = Adam's estimate gives F+(cx))-F|5'^^+F_(oo) = -0.0345 < 
in Eq. ( [7.15|) , which contradicts with the recent result from the lattice gauge theory. (See 



fig. 2.) There is a disagreement in the sign. 

(3) Numerical evaluation in the generalized Hartree-Fock approximation 

In this subsection we present various results obtained by numerical evaluation. The 
algorithm is simple. With given /iL, m//i, and 6es, we start with an initial fii/fi. Then Eq. 
( [7. 1| ) is solved numerically to find f{pw)- With this f{pw), a new is determined by 
([7.2|). We have a mapping 

/^i fiPw) /ii • (7.18) 

We repeat this process until the output fii coincides with the input /ii within required 
accuracy. With /ii being fixed, the condensate is evaluated by (|7^ ). 

In fig. 2 the m dependence of the condensate is plotted at = with various values of 
T/fi. The lattice result from ref. is also plotted for comparison. Our result agrees well 



with Tomachi and Fujita's evaluation by the Bogoliubov transformation. [|T6[ The agreement 
with the lattice result is modest. Note that the subtraction of condensates in free theory in 
each regularization scheme is crucial. 

There appears a singularity in the m dependence of the condensate when 6 is close to 
TT. In fig. 3 the condensates are plotted with various values of 6. The discontinuity appears 
at m/fi = .44 for 6 = n, and at m/ fi = .40 for 6 = .957r. 

The discontinuity persists so long as the temperature is lower than Tc ~ .12 fi. The 
condensate at various values of T is depicted in fig. 4. 

The origin of the discontinuity is understood as follows. 0] For fiL ^ 1 Eq. (|7.1| ) 
becomes 

^ ^ +^P^-^^^^^cos(2w-M|/(Pu^) = e/(M . (7.19) 



(27r)2rfp^ 4 "^"^ 47r2 
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Figure 2: The mass dependence of the chiral condensate in the = 1 model at = with various 
values of temperature T . In the figure T is in a unit of [i. The lattice data is by de Forcrand et 
al.[p^] The additional lattice data point at m/ ^ = 2 was provided by the authors of ref. [^^. The 
condensate for T/// < 0.1 is essentially the same as that at T = 0. The curve for T j [i = 0.003 is 
consistent with the result by Tomachi and Fujita.[16| 



The potential term dominates over the kinetic energy term. Suppose that ^efr = tt. 
Then the wave function is sharply localized around the absolute minimum of Vipw) = 
(nfipw)'^ + mfiie^ cos{2npw) ■ There is always a solution for which /x^ > 2e^mfii. In this 
case V{pw) is minimized at pw = and fii = v^/i^ + rd^e^^ — me"' . There is another 
solution for raj ^ > 0.435, or /i^ < 2e"'m^i. V{pw) is minimized at pw = Pw where 
2Trpw = (2m/iie'''//i^) sin27rpvy- The boson mass is given by /ii = y/JjP~+~aF — a where 
a = me"' cos2ttpw- The second solution has a lower energy density and corresponds to the 
vacuum. Hence at T = L^^ = 0, there appears a discontinuity at mc//i = 0.435. 

At finite temperature the critical value mc/ fi is determined numerically, mc/ = 0.437 
at T/fi = 0.03. The discontinuity disappears for T/fi > 0.12. As noted in ref. |jT2| there 
are two possible scenarios. In the full theory the discontinuity may persist, but with a 
universal m^/fi independent of T/yU. Or the discontinuity may be smoothed by higher order 
corrections. At the moment we do not know for sure which picture is right. 
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Figure 3: The mass dependence of the chiral condensate in the = 1 model near = tt at 
T/^ = 0.003. There appears a discontinuity above 9 ~ .95 tt. The condensate at ^ = is also 
displayed for comparison. 

8. Degenerate fermions 

When all fermion have degenerate masses {iria = m) and obey the same boundary 
conditions aa = a, the Laplacian and potential Vn{pw^ v) the eigenvalue equation 
( p.25|) become 



N-l Ql 



B 

Boson masses Hi and /i2 = 



N-l 



1 S ^^"^ 



dipb 



{liLf , NmLB ^ 



■Pw 



71 



cos (v^a - 



27rp 



w 



a=l 



N 



■ = Hn are determined by ( [4.5| ) - ( ^78| 
= fi^ + R , fil = R 



R=^^{cos [ipa — j )/ 



(8.1) 



When fiL, fiiL, fi2L ^ 1, the potential is approximated by 



12) 



l/N {N~l)/N 



' a=l 



N 



2ttpw 



(8.3) 
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Figure 4: The mass dependence of the chiral condensate in the = 1 model at = vr at various 
values of T/fi. For T//i < 0.12 the discontinuity remains, whereas the transition becomes smooth 
for r//i > 0.12. 



If fermion masses are small m/yU <^ 1, the first term in the potential dominates over the 
second. The pw dependence of the wave function is the same as in the massless case, and 
one can write f{pw, = e~'^'^^^'vv/'^^ f^Lp). With the aid of the truncation formula ( |6.2D the 
equation is reduced to 



{-A^ + «:F;vM}/M = e/(¥.) 



N 



N 



K 



tt{N 

The boson mass fi2 is given by 
2 SttttiB 



—-mLBe-^"'^'^ , ^^(y;) = - ^ cos<^„ 

^) a=l 



•4) 



l-K^{N-l) 



K {{ COS (fa )) f 



1.5) 



In the L ^ oo or T ^ limit, the wave function has a sharp peak at the location of the 
minimum of Voo{pw, ( p.3[ ) or Fn{^p) in ( |8.4| ). We examine the location of the minimum 
of the potential. 



(1) Potential 



(a) = 2 
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In the two flavor case (/^i = y?, = ^ — and the potential takes the form of 
Voo (pw, v) = -7-^ I {.T^l^Pwf - 4:6^ mfiy^fil^"^ cos (npw - ^) cos (v? - ^) I (8.6) 



Ait 



or 



F2(v9) = -2cos^cos(v9-^) . (8.7) 

The form of the potential Voo suggests that the anomalous behavior analogous to that in the 
= 1 case may develop near 6 = 71 with m = 0(/i). Note that at = ±7r the truncated 
equation ( |8.4| ) is not valid as the potential F2{ip) = 0. The pw degree of freedom must be 
retained. 

F2{ip) has a minimum at 



cn = - , = 9-271 

^ 2 ' 



9 + 71 
271 



where [x] denotes a maximum integer not exceeding x so that \9\ < tt. Notice that the 
location of the minimum discontinuously changes at 6' = vr {mod 27r). 

(b) = 3 

Let us examine the potential Fs{(f). 

F3{(p) = - cos{(pi + (p2 - 9) - cos (pi - cos (p2 ■ (8.9) 
The location of the minimum can be easily found analytically. There are six distinct 



stationary points: v^i = v^2 = \(^,\{(^ ± 27r), or {ipi,{p2) = {9 + 71, 9 + 71) 



3 ' 3 



71, 



{—9,9 + vr). The global minimum is located at fi = (p2 = The periodicity in 9 is 
271. The location of the minimum jumps from {Lpi,Lp2) = (+|7r, +|7r) to {—^71, —^71) at 
9 = n (mod 27r). The minimum is always located at \ipi\ = \ip2\ < ^tt. 

(c) General 

First notice that 

Fn{(Pu - ■ ■ ,(PN~1;6) = - XI COS - COS (6* - J2 ^cl) 

a=l a=l 
= FN-l{(pir ■ ■ ,^N-2;0 - (fN-l) - COSipN-l ■ (8.10) 



For = 3 we know the minimum is located at (fi = (p2 = ^0- 
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Consider = 4. We denote the minimum of F4 by (ai, 02, as). Fix the value of (ps and 
consider F4(97),^3fixed = G^i^ipi, ip2). Denote the minimum of by (61, 62) where bj = bji^ip^). 
It follows from ( |8.1U| ) that (&i,&2) is the minimum of F^{Lpi^ip2'-,d ~ fs)- The result in the 
= 3 case implies that 61 = 62- Since the minimum of H^^ip^) = G4^{bi[ip3,],b2[(p3]; fs), 
which we denote by c, is the minimum of F^lip), we have [01,02,03] = [61(c), 62(c), c]. In 
particular, oi = 02 as 61 = 62. We repeat the argument with the value of v?i kept fixed, to 
obtain 02 = 03. Hence Oi = 02 = 03, i.e. the minimum of F4^{ip) occurs at ip>a = 'P- 

By induction we conclude that FN{<f) has the minimum at = (a = 1, ■ ■ ■ , A^ — 1). It 
is easy to find the location of the minimum of F^lip) = Fn{ip>, ■ ■ ■ ,Lp;6). From the symmetry 
fN = 6 ~ {N ~ l)(p = (fi {mod 27t), or (f = 9/N {mod 2tt/N). Direct evaluation of F]^{if) 
shows that the minimum of Fiy{{p) is attained at = 9/N. 

(2) Boson masses and condensates 

In a few limiting cases boson masses and chiral condensates can be determined analyti- 



cally. In this subsection we suppose that fermion masses are small m -C /i and analyze (8^). 
The wave function f{<p) is determined by two parameters, 6 and k. The chiral condensates 
are related to the boson mass by 

{lpai'a)e = ii'ai'aYe - ii'alpaYfree 

i^ai^aYe = ■ (8.11) 



(a) = 2 

We suppose that Q ^ n. Eq. ( |8.4| ) becomes 
{ d^ 



dip" 



kq cos (v' - 2) f •^(V^) = ^ •^('^) 



^- 4mL rM/2D^,. rM/2„-7r/2uL „„„ ^ 



Ko = 2k cos \e = 5(/iiL)^/2^(/i2L)^/^e-"/2'^^ cos - . (8.12) 



It is easy to see 



COS ^6* = k(1 + cos6') for ^ 1 
cosy.))^=<; (8.13) 
cos ^9 for ^ 1 . 



27 



and accordingly [i2Ll\pl'K = kq or for kq -C 1 or kq S> 1, respectively. Hence one finds 
for m -C 



1^2 



' 4:^/2 m cos \e e-'^/2M^ for yuL < 1 
4v^m cos \e (^-i^) ^''^ for /iL > 1 > (m cos \Qfl^ii^l'^L 



2" \ 

\ 2/3 



(8.14) 



f 2e^ m;ui/2 cos \e j for (m cos \efl^[i^l^L > 1 



Coleman obtained /i2 oc m^/^yU^/^ cos^^^ |i9 in Minkowski space-time long time ago,]^ but the 
overall coefficient was not determined. Note that if the massless limit is taken with a fixed 
L, then /i2 = 0(m). 

At = and L — >■ oo, m dependence of /i2 has been determined in the lattice 
gauge theory. |p2| In this limit the formula ( ^.14|) leads to /i2/e = l^l^e'^^'^Tx^^l^imj e)^!'^ = 
2.163(m/e)^/^ where e and ce are the coupling constant and Euler's constant, respectively. 



Smilga showed that the exact coefficient should be 2.008 for m/e ^ The lattice simula- 

tion is done for m/e < .5. The data supports the m?^^ dependence and indicates a coefficient 
between the two numbers mentioned above. 



(b) iV > 3 flavor 



When K ^ 1 in Eq. ( ^.4| ), the wave function has a sharp peak at the global minimum 
= e/N, and therefore (( e'^-^ )) f = e'^^^ . 

For K <^ 1 we solve ( |8.4| ) in a power series in k. To 0(1) a plane wave solution 
u{(f;n) = e™^'''^"' i-^^iv-iviv-i gatisfles —A%u = eu where e is positive semi-deflnite and 
vanishes only if ni = ■ ■ ■ = riN^i = 0. Hence to 0(k°) f^^^ = 1. To flnd 0(/t) correction, we 
note A^F/v(v2) = —Fiy{{p), from which it follows that / = 1 — nF^^ip). Hence 

; for K -C 1 

cosy.J)^=<| 9 (8.15) 
' cos — for K 1. 

Notice that it is independent of 6* for k < 1. It follows that [N/87i'^{N - l)Y/^fi2L = k for 
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K <C 1 and = [kcos{6/N)]^^'^ for k 3> 1. The boson mass is determined from 

8N \i/2 



A^- 1 



m e"'^/^/^^ for /iL < 1 



N 



AT \ 1/2 /p7,, r \ 1/N 

j f j for /iL > 1 > m^/(^+i) L 



.16) 



n X Ar/{Ar+i) 

2e^ m/xV^ cos ^ J for m^/(^+i)/ii/(^+i)L > 1 



9. General fermion masses 



(1) = 2 

It is of interest to know the boson masses when the fermion masses are not degenerate. 
In the two flavor case an analytic expression is obtained for <^ in the L — oo limit. 
Start with (|3.25|) , or more conveniently the truncated equation (|6.3| ): 



rf2 _ 2L^, 



dip'' 



--/^A'^ COS J = eo/(¥') • 

a=l ) 



(9.1) 



Note that ipi = ip and (p2 = 9 — cp. In terms of Ra = SnmaBaL ^e'^^'^'^^i cosipa ))/, 

■2 ^ i?i + 



+ 



B{fi2L) 



±(i?i-R2)/2At2 



(9.2) 



In the L ^ oo limit, Ba oc L so that (( cosy9 ))/ = cos(/)min where the potential term in 
!T]) has a minimum at v^min- In general at the minimum of a function g{(p) = —a cosy? — 
psin^, e*^ = {a + il3)/^o?^TW- Hence, 

(irLaBaf + mim2-Bi-B2 cos 6 



R, = 87i^ 

y (mi5i)2 + (ms^a)^ + 2mim2BiB2 cos ^ 

where Ba = Ba/L. Solving (|9^ ) and (|9.3|) to the leading order in rria/fJ,, one finds 



(9.3) 



D _ 0^47/3,2/3 mg + mimacosg 

{mf + 1712 + 2mim2 cos c')^/'^ 



Consequently, 



/i2 = e'^"'^^jJ^^^{m'l + 777,2 + 27?li7772 COS ( 



a/3 



(9.4) 



(9.5) 
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Figure 5: The location of the global minimum of the potential (9.7) in the A'^ = 3 model with 
general fermion masses. As 9 changes from —it to vr, the minimum in the figure moves from left to 
right. The values for the fermion masses, rria's, are in a unit of fi. The location of the minimum 
discontinuously jumps at 9 = ibvr for (mi, m2, ms)//^ = (0.01,0.01,0.1) or (0.01,0.011,0.1), but 
makes a continuous loop for (0.01,0.02,0.1). 



and 

(^1^1 Ye 
In the symmetric case mi 



e47/3^2/3 



rrii + 1712 cos 9 



(9.6) 



27r(mi +1712 + 2mim2 cos 61)1/3 \^ 777,2 + mi cos^, 
1712 = ( |9-5|) reduces to ( ^.14|) . 
Observe that there is no singularity at 6^ = tt in a generic case mi 7^ 1712- The boson 
mass and chiral condensates are smooth functions of 6 with a period 27r. The singularity 
appears only when mi = 1712 7^ in the two flavor case. In a special case mi = but m2 7^ 



(m2 < /i), ((^1^1 (^/'2^2 )0)=-(27r) 



'^{e^"' fi^7n2Y^^{cos9, 1) in this approximation. 



(2) = 3 

The three flavor case mimics physics of four-dimensional QCD. Three fermions, which 
one may call "up", "down", and "strange" quarks, have different masses; mi ~ m2 <^ m^. 
We would like to see how the asymmetry in masses affects chiral condensates and boson 
masses. 

We concentrate on the L ^ 00 (T = 0) limit in the truncated theory. As in the N = 2 
case one needs to find the location of the minimum of the potential 



TT 



9{^l,^2;0) = lim ITJ^^Viif) = - J^ITT'aB a cos Lf a 



(9.7) 



a=l 
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where (ps = 9 — ipi — ip2- Here B^s are determined by the boson masses /Xq's and eigenvectors: 

a=l 

Ra = 87rm,5,cos(/^r" • (9i 




is the location of the minimum of g{(f). The set of equations ( |9.71) and ( p.8|) must be 
solved simultaneously. Chiral condensates are given by 

( ^a^a = -2Ba COS </.r + — • (9.9) 

71 

When fermion masses are degenerate, v^™" = ^0, as shown in Section 9. At = tt 
the location of the minimum changes discontinuously, which induces singular behavior in 
physical quantities. We show that the singularity disappears if the asymmetry in fermion 
masses is sufficiently large. In ref. |11| this problem was analyzed by examining gif), but 



without solving ( |9.8|) . 

When mi = m2 < rris, ifii = if2 at the minimum of g{(f). At 6 = 0, (pi = ip2 = 0. As 
6 increases, = f2 also increases. At = vr it reaches <fi = if2 = fc whose value depends 
on m^'s. As gets bigger and bigger with mi = m2 kept fixed, (pc approaches ^vr. For 
instance <pc = (0.467,0.486,0.499)77 for nii/fi = 0.01 and ms/fi = (0.02,0.03,0.1). As 6 
exceeds vr, the minimum jumps to ipi = {p2 = —(fc and returns to ipi = ip2 = at 9 = 27r. 
The singular behavior at 6* = vr remains. This is expected as ^ mi = 1712 corresponds to 
the two fiavor case. 

Now we add a small asymmetry in the light fermions. Fig. 5 depicts the location of 
the minimum of g{ip) as 6 varies from — vr to vr. A small asymmetry in mi and 1712, does 
not change the behavior near ^ = 0, but significantly affects the behavior near 6 = ivr. At 
{1711,1712,7713)/ fi = (0.01,0.011,0.1) the minimum at 6* = tt is very close, but not quite equal, 
to (7r,0). At (1711,7112,7173)/ fi = (0.01,0.02,0.1) the minimum at 6' = vr is located at (vr, 0), 
and the singularity in physical quantities at = vr disappears. 

In fig. 1^ we have plotted boson masses /ii, fi2, and fi^ as functions of 9 with given fermion 
masses. They correspond to m^', m^, m.,^ in QCD. For mi = 7712 = 7713 <^ yU, /ii ^ /i2 = /is. 
When mi = m2 < 7713 <^ fi, all fi^s are different. The 6 dependence of each /i^ has similar 
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Figure 6: Boson masses in the N = 3 model at T = with degenerate fermions mi = m2 = 
and with nii = m2 < m^. In the former case 1^2 = Fermion masses are in a unit of fi. 



behavior, fii/fi, /i2//^, and fis/^ vary by 0.003, 0.007, and 0.1 in magnitude. The mass of 
the hghtest boson, m^, has the most 6 dependence. 

In fig. ^ the 9 dependence of the mass of the hghtest boson is plotted for various values 
of fermion masses. The cusp at = vr persists so long as mi = 1712, but a small asymmetry 
in mi and m2 changes it to smooth dependence. The mass at 6* = increases as the third 
fermion gets heavier as expected, whereas it decreases at = vr. 

The 6 dependence of the location of the minimum of g{^p) and the value of the mass 
/i3 of the lightest boson induces nontrivial 6 dependence in the chiral condensates {ip a'^a )e- 
In fig. 1^ a chiral condensates for (mi, m2, ms)//! = (0.01,0.01,0.01) and (0.01,0.01,0.1) are 
depicted. In both cases there appear cusps aX 9 = -k. Notice that the magnitude of the 
condensates at ^ = is insensitive to fermion masses. This is true only if the "free theory 
background", {ipi^Yfree is subtracted in the definition of the condensate ( p.9| ). 

There appears, however, a big difference in the 6 dependence of the condensates. When 
(mi, m2, m3)/yU = (0.01,0.01,0.1), the third fermion ip^ is much heavier than the other two. 
The condensate ( V'sV's )e is more or less independent of 9, which is expected as the vacuum 
structure is mainly determined by light fermions. 

When a small asymmetry in light fermions is added, condensates suffer a big change. See 
fig. p. With (mi,m2,m3)/yU = (0.01,0.011,0.1) the 9 dependence in {ipiipi)0 is enhanced, 
whereas {ip 2^^2)0 develops a dip near 9 = tt. A small asymmetry in mi and m2 induces a 
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Figure 7: The mass, /X3, of the hghtest boson in the A'^ = 3 model at T = for various values 
of fermion masses. Fermion masses are in a unit of /i. The cusp at ^ = vr disappears as a small 
asymmetry in the two light fermions is added. 

big difference in {ipiipi)g and {ip 2^^2)9 ^^ar 9 = n. The nonlinearity in Eq. and (|9.8| ) 
gives rise to such sensitive dependence. 



It is interesting to recognize the similarity between the potential g{(p) in (|9.71) and the 
effective chiral Lagrangian proposed by Witten to describe low energy behavior of four- 
dimensional QCD.[^ In Witten's approach 

^Witten^j;) = /2| - -Tr M{U + U^) + -^{-i\ndetU -ef] (9.10) 

where U is the pseudo-scalar field matrix and M = diag{my_,mii,ms) is the quark mass 
matrix. The second term represents contributions from instantons. k is 0(1) in the large Nc 
(color) limit. 

The fact that m^, ^ m1,m'j^,m'^ implies that k/Nc^ ma- Diagonalize U and denote 
it by (im(yf(e*'^^, e**^^, e*"^^). As the second term dominates over the first, l]0a = ^ to the 
first approximation. With 03 eliminated, ^^'**'^°([/) takes the same form as g{^p) in (|9.7|) . 
Consequently both models show qualitatively similar behavior. Indeed, Witten has argued 
that a small asymmetry in m„ and rrid, in addition to large asymmetry rris ^ mu, rud removes 
the singularity of physical quantities in 6 at vr, which is exactly what we are observing in the 
= 3 Schwinger model. However, it should be noted that the coefficients in the potential 
g{^p) have extra fermion mass dependence coming from the factors Ba, which have significant 
effects near 6* = vr. 
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Figure 8: The 6 dependence of the chiral condensates in the N = 3 model at T = 0. Two cases, 
{mi, 1712, ms)/iJ, = (.01, .01, .01) and (.01, .01, .1), are displayed. In the latter case the condensate of 
the heavy fermion has little dependence on 6, whereas the condensates of the light fermions show 
large dependence. 

10. Summary 



In this paper the massive N-flavor Schwinger model was analyzed in the generalized 
Hartree-Fock approximation. Dynamics of the zero-modes is determined by the Schrodinger 
equation for N degrees of freedom. The potential term in the Schrodinger equation depends 
on the boson spectrum in the oscillatory modes. The boson mass spectrum in turn depends 
on the ground state in the zero-mode sector. The ground state of the two sectors must be 
determined self-consistently. 

We have evaluated the boson spectrum and chiral condensates in the iV = 1, 2, 3 
models. In the A'^ = 1 model we have found anomalous dependence of physical quantities on 
the fermion mass near ^ = tt at low temperature. In the N — 3 model physics near 9 — tt 
is very sensitive to the small asymmetry in fermion masses. Chiral dynamics in the N — 3 
model resembles with that in QCD in four dimensions. 
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Figure 9: The 9 dependence of the chiral condensates in the = 3 model at T = 0. Small 
asymmetry is added to masses of light fermions. The effect is minor near 6 = 0, but the condensates 
of the light fermions near 9 = tt are significantly affected. 



A. Correlation functions 

Green's function for a scalar field with a mass n on S^, excluding the zero mode, is 



In the massless case (/x = 0) it is given by 



cos 



2mnx 



(A.l) 



A(x;0,L) = --^-ln2fl- 



2tix 



4:71 

L 



cos 



2Trx 



for 



27rx 



L 



< 1 



In terms of 



I- 7 n v-^ COS 2'rinb 
I[s;a,b = 2^ 



(A.2) 



(A.3) 



n=l 



2TrA{x;fx,L) = J[ |; /iL/27r, x/L ] and 5(/iL) = exp { - /[ i; /xL/27r, 0] + /[ i; 0, 0]}. For 
161 < 1 



I[s;a,b] = / 
Jo 



cos 2nbt 1 



+ 2 sin STT / (it 

J a 



cosh 27r6t 



(t2-a2)^(e2^*-l) 



r(.) 



i(27r|a6|)- 



a 



2a2^ 
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+2 sin sn dt- 

J a 



cosh 27r6t 



(t2-a2)^(e2-*-l) ■ 

The above formula is vahd for an arbitrary s by analytic continuation. 
It then follows that 

cosh /iXM 



(A.4) 



27rA(x; /i, L) = Ko{\ij,x\) h 2 / du 

fiL Ji 



(A.5) 



Recalling that Ko(z) ~ — In 2 — 7 + In 2 for z -C 1 and ~ y ^ e ^ for z ^ 1, and noticing 
for /iL > 1 

cosh /iXM r 



2 / c/m 



TT 



-fi(L—x) 



1 (e'^^" - - 1 V 2/i(L 



+ 



TT 



X] 



2/x(L + x) 



-/^(L+a::) 



(A.6) 



one finds 



g27rA{a;;/i,L) ^ . 



X 



1 for /iL > 1 , - < 1 , /ix > 1 



2e~'^ x 

for uL > 1 , — <^ 1 , ax . 

fix L 



(A.7) 



B. Normal ordering and Bogoliubov transformation 

In the subsequent discussions, we make frequent use of identities: (1) : := e^^, 
(2) = e^I^'^l e^+^ = el^'^l e^, and (3) : : : := el^^'^^l : e^+^ : where A+ and 

y4~ denote the annihilation and creation operator parts of A, respectively. With massless 
fields (13) 

4:71 

We also note that 

1 ( p+iTTx/L —i-Kx/L ^ 

- ) - \ - I = e*''^'/^ 6r (x) 

1 _ Q+27,i{x+ie)/L - ^ Va; + ze ^ yZ ^ 360 + ' ' ' J ' *^ ' 

We have seen that boson fields become massive due to the Coulomb interaction and 
fermion masses. When this happens, the vacuum also changes and in two dimensions non- 
vanishing chiral condensates result. 
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It is most convenient to work in the Schrodinger picture. A boson field is generically 
denoted by with its conjugate n(x). On a circle they are expanded as 

1 



E 



n(x) = E \/^{cn(/^)e^^"^ - ct (/.)e-^-} 



(B.3) 



where Pn = 2im/L and ^^(yu) = y'p^ + yU^. Annihilation and creation operators c„(/i) and 
c.|^(yu) are defined with respect to a mass /i. The left- and right-moving modes, (p+^x) and 
of = + 0_ are defined by the n < and n > components in ( [B.3| ). In the 
massless limit they corresponds to (|2.2|) in the Schrodinger picture. 
Sometimes we need to treat 0±'s separately. One finds 



(t)±{x) = l(f){x) ± -5- 

/ ZL Jo 



dyF{x -y;fi) U^{y) 



i ( 

^(^;/^) = E--7-Y^ 



— 2ninx/L +2-Kinx/L 



} 



(B.4) 



1l ^nifJ-) 

Note that the definition of (f)±{x) depends on the reference mass fi as opposed to that of 0(x 
and Il^{x). In the massless theory F'{x; 0) = L5l{x) — 1 so that n(x) = 0+(a;) — 0'_(x). 
We have an identity among c„(/i)'s with different /x's: 

cosh6'„(/ii;/i2)c„(/i2) + sinh6'„(/ii;/i2)cL„(/i2) 



C„(/il 



"cosh 6'„(/ii;/i2) 
.sinh 6'n(/ii;/i2) 



t^n(/il 



± 



\ t^n(/i2) \ t^n(/il) 



u;„(/i2 



(B.5) 



In other words the change in the boson mass induces a Bogoliubov transformation. The 
vacuum with respect to a boson mass is defined by Cn(/i) |vac; fi) = 0. 

In our formalism fermion fields are first bosonized in the interaction picture defined by 
massless bosons. Boson fields then acquire masses and the vacuum is redefined. In particular 
boson fields are normal-ordered with respect to physical boson masses. One useful relation 
is 



a 



exp 



2 oo 



(B.6) 



from which it follows P, ^ 
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B{i,L) = exp - 5: (- 



^2 n 



= ^ exp <^ 7 + — - 2 / — — r \ . B.7 

In the following we make use of simplified notation: ^±{x) refers to the massless field, 
— <^n(0) and Qn — ^n(0;/i). iV^[A(c, c''")] denotes that the operator A is normal-ordered 
with respect to c„(/u) and c|^(/i). Some useful identities are 

A^o[ 6*°'^^'^^+*^'^^^^ ] /tV^J gja</'W+i/3</'(y) j _ g-a/3{A(a;-2/;/<,L)-A(i-2/;0,L)} (B.8) 



= exp s — X! — r 2q;+q;_ cosh ^„ sinh ^„ + (a^ + a^) sinh^ Qn \ 

= 5(/iL)("^+'^^)/^'' g-a/3{A(x-j/;M,I')-A(j;-j/;0,L)} g-|(a2+/32)/i(0;M,L)-a/3?i(j;-j/;M,L) 



oo 



C. Useful identities for zero modes 

(1) Free massless fermion 

In a theory of a free massless fermion the Hamiltonian in the bosonization method is 

given by 

H = ^(V + / dx^i^' + (C.l) 
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where = 0+ + 0_ represents a massless boson defined in ( p.2D . The vacuum state is 
1^^^^ ) = |72 = ), i.e. p l^vac ) = 0. Note that ( e'°"^ ) = 6^,0- Since q{t) =q + Airpt/L, 



' ^iaq(t) ^ibq{0) ^ ^^^^ ^-2nia'^t/ L 



^y[g^ag(^)g^65(0)]^ = ^,+5 e-^^-'W/i . (C.2) 

(2) QED2 

If all fermions are massless, the model is solved in the operator form in the Heisenberg 
picture. The solution to ( p.4|) in the iV-fiavor model is 



^wytj = Bv^cos/xt + -^iVsm/it 
Pw{t) = Pw cos fit -Q'yysinut 

TCfiL 

Airt 1 

Pa{t) = Pa 

qa{t) = qa + -^Pa 

Pa{t) = Pa . (C.3) 

Here Q'^^ and q'^ are defined in (|3.3| ). 
In the = 1 model 

^ ^ieq+il3{ew+2TTp)+i-ypw ^ ^ ^iee ^-TT/iLlS^ /4-{-y-2TTef /4ttiiL 

It follows from (^) and (^) that 

( e-^(^)e-^'^(°) ) = exp L{a + (3)9 - | - '^e'A . (C.5) 
I fiL fiL J 
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